Gyrokinetic Landau collision operator in conservative form

A gyrokinetic linearized exact (not model) Landau collision operator is derived by transforming the symmetric and conservative Landau form. The formulation obtains the velocity-space flux density and preserves the operator's conservative form as the divergence of this flux density. The operator contains both test-particle and field-particle contributions, and finite Larmor radius effects are evaluated in either Bessel function series or gyrophase integrals. While equivalent to the gyrokinetic Fokker–Planck form with Rosenbluth potentials [B. Li and D. R. Ernst, Phys. Rev. Lett. 106, 195002 (2011)10.1103/PhysRevLett.106.195002], the gyrokinetic conservative Landau form explicitly preserves the symmetry between test-particle and field-particle contributions, which underlies the conservation laws and the H theorem, and enables discretization with a finite-volume or spectral method to preserve the conservation properties numerically, independent of resolution. The form of the exact linearized field-particle terms differs from those of widely used model operators. We show the finite Larmor radius corrections to the field-particle terms in the exact linearized operator involve Bessel functions of all orders, while present model field-particle terms involve only the first two Bessel functions. This new symmetric and conservative formulation enables the gyrokinetic exact linearized Landau operator to be implemented in gyrokinetic turbulence codes for comparison with present model operators using similar numerical methods.United States. Department of Energy (Contract DE-FC02-08ER54966

Momentum transport, radial electric field, and ion thermal energy confinement in very high temperature plasmas

Thesis (Ph. D.)--Massachusetts Institute of Technology, Dept. of Physics, 1998.Includes bibliographical references (p. 213-238).Darin R. Ernst.Ph.D

Stationary density profiles in the Alcator C-mod tokamak

In the absence of an internal particle source, plasma turbulence will impose an intrinsic relationship between an inwards pinch and an outwards diffusion resulting in a stationary density profile. The Alcator C-mod tokamak utilizes RF heating and current drive so that fueling only occurs in the vicinity of the separatrix. Discharges that transition from L-mode to I-mode are seen to maintain a self-similar stationary density profile as measured by Thomson scattering. For discharges with negative magnetic shear, an observed rise of the safety factor in the vicinity of the magnetic axis appears to be accompanied by a decrease of electron density, qualitatively consistent with the theoretical expectations. © 2012 American Institute of Physics.United States. Department of Energy. Office of Fusion Energy Science

Characterization of density fluctuations during the search for an I-mode regime on the DIII-D tokamak

The I-mode regime, routinely observed on the Alcator C-Mod tokamak, is characterized by an edge energy transport barrier without an accompanying particle barrier and with broadband instabilities, known as weakly coherent modes (WCM), believed to regulate particle transport at the edge. Recent experiments on the DIII-D tokamak exhibit I-mode characteristics in various physical quantities. These DIII-D plasmas evolve over long periods, lasting several energy confinement times, during which the edge electron temperature slowly evolves towards an H-mode-like profile, while maintaining a typical L-mode edge density profile. During these periods, referred to as I-mode phases, the radial electric field at the edge also gradually reaches values typically observed in H-mode. Density fluctuations measured with the phase contrast imaging diagnostic during I-mode phases exhibit three features typically observed in H-mode on DIII-D, although they develop progressively with time and without a sharp transition: the intensity of the fluctuations is reduced; the frequency spectrum is broadened and becomes non-monotonic; two dimensional space-time spectra appear to approach those in H-mode, showing phase velocities of density fluctuations at the edge increasing to about 10 km s−1. However, in DIII-D there is no clear evidence of the WCM. Preliminary linear gyro-kinetic simulations are performed in the pedestal region with the GS2 code and its recently upgraded model collision operator that conserves particles, energy and momentum. The increased bootstrap current and flow shear generated by the temperature pedestal are shown to decrease growth rates, thus possibly generating a feedback mechanism that progressively stabilizes fluctuations.United States. Department of Energy. Office of Fusion Energy Sciences (Award DE-FG02- 94ER54235)United States. Department of Energy. Office of Fusion Energy Sciences (Award DE-FG02-94ER54084)United States. Department of Energy. Office of Fusion Energy Sciences (Award DE-FG02-08ER54984)United States. Department of Energy. Office of Fusion Energy Sciences (Award DE-FC02-04ER54698

The effects of main-ion dilution on turbulence in low q95 C-Mod ohmic plasmas, and comparisons with nonlinear GYRO

Recent experiments on C-mod seeding nitrogen into ohmic plasmas with [subscript q]95 = 3.4 found that the seeding greatly reduced long-wavelength (ITG-scale) turbulence. The long-wavelength turbulence that was reduced by the nitrogen seeding was localized to the region of r/a≈0.85, where the turbulence is well above marginal stability (as evidenced by Q[subscript i]/Q[subscript GB]≫1). The nonlinear gyrokinetic code GYRO was used to simulate the expected turbulence in these plasmas, and the simulated turbulent density fluctuations and turbulent energy fluxes quantitatively agreed with the experimental measurements both before and after the nitrogen seeding. Unexpectedly, the intrinsic rotation of the plasma was also found to be affected by the nitrogen seeding, in a manner apparently unrelated to a change in the electron-ion collisionality that was proposed by other experiments.United States. Dept. of Energy. Office of Fusion Energy Sciences (Award E-FG02-94-ER54235

Poloidal asymmetries in edge transport barriers

Measurements of impurities in Alcator C-Mod indicate that in the pedestal region, significant poloidal asymmetries can exist in the impurity density, ion temperature, and main ion density. In light of the observation that ion temperature and electrostatic potential are not constant on a flux surface [Theiler et al., Nucl. Fusion 54, 083017 (2014)], a technique based on total pressure conservation to align profiles measured at separate poloidal locations is presented and applied. Gyrokinetic neoclassical simulations with XGCa support the observed large poloidal variations in ion temperature and density, and that the total pressure is approximately constant on a flux surface. With the updated alignment technique, the observed in-out asymmetry in impurity density is reduced from previous publishing [Churchill et al., Nucl. Fusion 53, 122002 (2013)], but remains substantial (nz,H/nz,L∼6). Candidate asymmetry drivers are explored, showing that neither non-uniform impurity sources nor localized fluctuation-driven transport are able to explain satisfactorily the impurity density asymmetry. Since impurity density asymmetries are only present in plasmas with strong electron density gradients, and radial transport timescales become comparable to parallel transport timescales in the pedestal region, it is suggested that global transport effects relating to the strong electron density gradients in the pedestal are the main driver for the pedestal in-out impurity density asymmetry.United States. Department of Energy (DE-FC02-99ER54512)United States. Department of Energy (DE-FG02-06ER54845)United States. Department of Energy (DE-FG02-86ER53223)United States. Department of Energy (DE-AC02-09CH11466

Local and global Fokker-Planck neoclassical calculations showing flow and bootstrap current modification in a pedestal

In transport barriers, particularly H-mode edge pedestals, radial scale lengths can become comparable to the ion orbit width, causing neoclassical physics to become radially nonlocal. In this work, the resulting changes to neoclassical flow and current are examined both analytically and numerically. Steep density gradients are considered, with scale lengths comparable to the poloidal ion gyroradius, together with strong radial electric fields sufficient to electrostatically confine the ions. Attention is restricted to relatively weak ion temperature gradients (but permitting arbitrary electron temperature gradients), since in this limit a delta-f (small departures from a Maxwellian distribution) rather than full-f approach is justified. This assumption is in fact consistent with measured inter-ELM H-Mode edge pedestal density and ion temperature profiles in many present experiments, and is expected to be increasingly valid in future lower collisionality experiments. In the numerical analysis, the distribution function and Rosenbluth potentials are solved for simultaneously, allowing use of the exact field term in the linearized Fokker-Planck collision operator. In the pedestal, the parallel and poloidal flows are found to deviate strongly from the best available conventional neoclassical prediction, with large poloidal variation of a different form than in the local theory. These predicted effects may be observable experimentally. In the local limit, the Sauter bootstrap current formulae appear accurate at low collisionality, but they can overestimate the bootstrap current near the plateau regime. In the pedestal ordering, ion contributions to the bootstrap and Pfirsch-Schluter currents are also modified

Application of ECH to the Study of Transport in ITER Baseline Scenario-like Discharges in DIII-D

Recent DIII-D experiments in the ITER Baseline Scenario (IBS) have shown strong increases in fluctuations and correlated reduction of confinement associated with entering the electron-heating-dominated regime with strong electron cyclotron heating (ECH). The addition of 3.2 MW of 110 GHz EC power deposited at ρ~0.42 to IBS discharges with ~3 MW of neutral beam injection causes large increases in low-k and medium-k turbulent density fluctuations observed with Doppler backscatter (DBS), beam emission spectroscopy (BES) and phase-contrast imaging (PCI) diagnostics, correlated with decreases in the energy, particle, and momentum confinement times. Power balance calculations show the electron heat diffusivity χ[subscript e] increases significantly in the mid-radius region 0.4<ρ<0.8, which is roughly the same region where the DBS and BES diagnostics show the increases in turbulent density fluctuations. Confinement of angular momentum is also reduced during ECH. Studies with the TGYRO transport solver show that the model of turbulent transport embodied in the TGLF code quantitatively reproduces the measured transport in both the neutral beam (NB)-only and in the NB plus EC cases. A simple model of the decrease in toroidal rotation with EC power is set forth, which exhibits a bifurcation in the rotational state of the discharge.United States. Dept. of Energy (DE-FC02-04ER54698)United States. Dept. of Energy (DE-FC02-08ER54966)United States. Dept. of Energy (DE-AC03-09CH11466)United States. Dept. of Energy (DE-FG02-04ER54235)United States. Dept. of Energy (DE-FG0289ER53296)United States. Dept. of Energy (DE-FG02-08ER54999)United States. Dept. of Energy (DE-FG02-08ER54984)United States. Dept. of Energy (DE-FG02-04ER54461